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Abstract 

We report first-principles calculations of current-induced forces in Si atomic wires as a function of bias and wire length. We 
find that these forces are strongly nonlinear as a function of bias due to the competition between the force originating from 
the scattering states and the force due to bound states. We also find that the average force in the wire is larger the shorter the 
wire, suggesting that atomic wires are more difficult to break under current flow with increasing length. The last finding is in 
agreement with recent experimental data. 
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Current-induced atomic motion (electromigration) has 
always been of concern in microelectronics since it was 
found to be a major failor mechanism in aluminum con- 
ductors. ^ Such an effect is more pronounced the smaller 
the material dimensions as it is the case, for instance, 
for metallization lines in conventional circuits. Q It does 
not come as a surprise then that this effect has recently 
attracted considerable interest in nanoscience. In par- 
ticular, current-induced forces in atomic 1^, ^, ||, y, 0] 
and molecular wires ^ have been investigated. These 
systems show interesting physical properties which make 
them promising candidates for future electronic applica- 
tions. ||lo[| It was found that current-induced forces in a 
nanoscale wire can distort its local atomic structure and 
affect its stability according to the external bias and the 
wire length. |Q, ||, ||, |^, ||, Q Longer atomic wires, for 
instance, have been found to be more stable with respect 
to current-induced forces than shorter wires, 

Over the years considerable theoretical and experimen- 
tal work has been devoted to understand the microscopic 
origin of electromigration. ||l3| It has become clear, for 
instance, that this phenomenon depends strongly on the 
microscopic details of the self-consistent electric field that 
is created upon scattering of the electrons across the re- 
gion of interest. Despite the progress in the fundamen- 
tal understanding of current-induced forces, however, it 
is still unclear what determines, e.g., the resistance of 
atomic and molecular wires to electromigration effects or 
what is the dependence of these forces on the external 
bias and wire length. 

In this letter we explore the above issues using first- 
principles approaches. In particular, we study current- 
induced forces as a function of bias in Si atomic wires of 
varying length (we consider wire lengths of up to four Si 
atoms). We find that current-induced forces are strongly 
nonlinear as a function of bias due to the competition be- 
tween the force originating from the continuum of states 
and the force due to states in the discrete part of the 
spectrum. The force originating from the bare electrodes 
(we will call this force direct [M) varies almost linearly 
as a function of bias. We also find that the average force 
in the wire is larger the shorter the wire, suggesting that 
atomic wires of increasing length are more difficult to 
break under current fiow. The last finding is in agree- 
ment with recent experimental data. |I2] 



A schematic of the system investigated is depicted 
in the inset of Fig. |l|. It consists of a Si atomic wire 
sandwiched between two gold surfaces that we model 
with ideal metals (jellium model), jl^, |l^ The inte- 
rior electron density of the electrodes is taken equal 
to the value for metallic gold (r^ w 3). The elec- 
tric current is calculated using the method described in 
Refs. HF 
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The electron wave functions are computed 
by solving the Lippman-Schwinger equation iteratively 
to self-consistency in steady state. Exchange and cor- 
relation are included in the density-functional formal- 
ism within the local-density approximation. ||l^, |l^ The 
current is computed from the wave functions lip) of the 
electrode-molecule system. The force F acting on a given 
atom at position R due to the electron distribution as 
modified by the external bias is given by the Hellmann- 
Feynman-type of theorem developed in Ref. 0: il^ 
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The sum and integral in Eq. (Q) include spin variables 
also. The first term on the RHS of Eq. (|l|) is the usual 
Hellmann-Feynman contribution to the force due to lo- 
calized electronic states \tpi)- The second term is the con- 
tribution to the force due to the continuum of states. |^ 
It is calculated by constructing, for each energy in the 
continuum, square-integrable wavefunctions {iPa) in an 
energy region A 
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where ^ is a normalization constant and the ip^s are 
single-particle wavefunctions in the continuum, solutions 
of the Lippmann-Schwinger equation. ||^ The continuum 
integration a covers the part of the spectrum occupied 
by the electrons at a given bias. Finally, the total force 
on the atom includes a trivial ion-ion interaction. 

We start the calculations by first relaxing the atomic 
positions at zero bias. For all different wire lengths the 
relaxed Si-jellium surface bond length is about 2 a.u., 
and the relaxed Si-Si bond distance is about 4.2 a.u. In 
Fig. |l|(a),(b),(c) we plot the total force as a function of 
bias on each Si atom for wires composed of two, three 



and four Si atoms, respectively. The atomic positions 
have been fixed at the equiUbrium position at zero bias, 
and the atoms are labeled with increasing number start- 
ing from the closer to the left electrode (see Fig. 0). In 
Fig. |^(d) we plot the average force (sum of the forces on 
each atom divided by the number of atoms in the wire) 
as a function of bias and the force for a wire composed of 
a single Si atom. Positive force pushes the atom to the 
right, i.e., opposite to electron flow (the left electrode is 
positively biased). It is clear from Fig. [l](a),(b),(c) that 
the force on atoms is a nonlinear function of the bias 
(exept for the single-atom wire, see below and Rcf. |l|) 
while the average force Fig. ^d) saturates with the num- 
ber of atoms in the wire. We first discuss the nonlinear 
behavior of the forces and then discuss the average force. 

For a small external bias (0.01 V, i.e., linear-response 
regime), the current-induced forces in the 3-Si wire sat- 
isfy the zero-force sum rule while for the 2-Si and 4-Si 
wire all atoms are pushed opposite to the current flow 
with almost equal force as it can be inferred from symme- 
try considerations . On the other hand, strong nonlin- 
earities in the current-induced forces appear at voltages 
above 0.1 V. For instance, in the 2-Si wire the Si atom 
closer to the left electrode (Sii in the inset of Fig. |l|) is 
pushed against the electron flow for biases less than 0.5 
V and moves along the electron flow for larger biases. On 
the other hand, the second Si atom (Si2 in the inset of 
Fig. |l|) is pushed against the electron flow at all external 
voltages. Similar effects, involving different atoms in the 
wire, occur at even smaller voltages for the 4-Si wire. In 
the same vein, for the 3-Si wire the zero-force sum rule 
is not satisfied already at 0.1 V. It is interesting to note 
that for large biases the largest force occurs on the second 
Si atom from the left. The break-up of the wires is thus 
likely to nucleate from the bonds of this atom. In general, 
each atom experiences a force due to current flow that is 
nonlinear in the external bias. Such nonlinearities affect 
the atomic redistribution in the wire and eventually its 
resistance to current-induced rupture. In order to un- 
derstand this nonlinear behavior, we study the different 
contributions to the forces as a function of bias: (a) the 
contribution from the electrodes without the atoms in 
between (direct force |l^ ) , (b) the contribution from the 
continuum part of the spectrum, and (c) the part of the 
force originating from the discrete spectrum. The force 
originating from the ion-ion interaction does not depend 
on the bias and we thus not discuss it here. For the sake 
of simplicity we discuss the 3-Si wire case. Similar con- 
siderations are valid for other wire lengths. We plot in 
Fig. H the three different contributions considered. As 
expected, the direct force is almost linear with the bias. 
Small deviations from linearity appear due to the small 
deviations from linear decay of the electrostatic poten- 
tial close to the electrode surfaces. (|l5[ ^ Furthermore, 
due to this deviation, the force on the central atom in 
the wire is larger than the force on the two atoms close 
to the electrode surfaces at any bias (see Fig. ^). The 
latter atoms also experience a force of similar magnitude 



with deviations occurring at large voltages. From Fig. g 
is also evident that for biases larger than 0.5 V the major 
contribution to the force on each atom comes from the 
direct force. The nonlinearities in the total force then 
orginate from the competition between the force due the 
states in the continuum (scattering states) and the force 
due to the bound states (in this case, those states be- 
low the band bottom of the left electrode) . The force on 
the central atom of the wire (Si2) due to the scattering 
states pushes the atom to the left, i.e. along with the 
electron flow. This force is almost linear with the bias. 
The force on the same atom due to the discrete spec- 
trum, on the other hand is almost zero even at very large 
voltages. This behavior can be understood by looking at 
the extra charge localized on the atom at any given bias. 
This quantity can be estimated by integrating the charge 
around a given atom in a box centered on that atom with 
one face parallel to the electrode surfaces and whose size 
in the perpendicular direction is equal to the bond dis- 
tance between the atoms. The difference between this 
charge at zero bias and the corresponding charge at any 
given bias is plotted in Fig. |3| for the three different Si 
atoms. Apart from small fluctuations as a function of 
bias, the extra charge on the central atom is zero. The 
charge on this atom is thus practically constant at any 
bias. Since, in this system, the number of bound states 
does not change with bias, their contribution to the force 
is constant [see Fig. ||(c)]. On the other hand, the force 
due to the continuum [second term on the RHS of Eq. 
(1)], increases linearly with bias due to the larger energy 
integration in Eq. (1) for larger voltages. [ p8| | 

The continuum and discrete contribution to the force 
is nonlinear for the atoms closer to the electrodes (Sii 
and Sis of Fig. ||). For the Sii atom, for instance, the 
continuum force is almost constant up to about 1 V, 
while the force from the discrete spectrum increases in 
magnitude with bias in the same voltage range (the sign 
of the force corresponds to the atom pushed along with 
the current flow). For biases above 1 V, the continuum 
force increases in magnitude and the force due to bound 
states is almost constant. An opposite trend is observed 
for the Sis atom (see Fig. ||). The opposite trend can be 
explained again by looking at the extra charge on these 
atoms (Fig. ^). This charge is of similar magnitude but 
of opposite sign indicating that there is a charge trans- 
fer from the right electrode to the left electrode when 
current flows with consequent creation of a local electric 
dipole. |l|| 

While it is difficult to extract general trends in the 
overall resistance of the wires to electromigration by look- 
ing at the forces on each atom, the average force follows 
quite a simple trend as a function of bias and wire length. 
In Fig. |l|(d) we plot such a quantity for wires of two, three 
and four Si atoms. For comparison we also plot the force 
for a single Si atom. The average force is almost linear as 
a function of bias even for large biases (deviations occur 
above 2.5 V for the 2-Si wire). Furthermore, the force is 
larger the smaller the wire length and, in particular, it 
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is almost equal for the 3- and 4-Si wires. |^ This trend 
indicates that short wires are easier to break than longer 
wires, and the average force reaches a "bulk" value at 
wire lengths of only three atoms, i.e., wire lenghts of less 
than lOA. Similar trends have been observed in experi- 
ments with Au atomic wires. [|l2| 

We conclude this paper by studying the effects 
of current-induced atomic relaxations on the current- 
voltage (I-V) characteristics of atomic wires. It was found 
by Di Ventra et al. ^, |ll| that current-induced atomic 
relaxations in molecular wires do not substantially af- 
fect the absolute value of the current for large voltages 
and current densities. We find that this trend is also 
valid in the present case of atomic wires, suggesting it 
is a general trend valid for both atomic and molecular 
wires when current flow is mainly coherent. This is illus- 
trated in Fig. H where we plot the I-V characteristics for 
the 3-Si wire with and without current-induced atomic 
relaxations. In the inset of the same figure we show the 
relaxed atomic positions for each atom at selected biases. 
It is clear from Fig. ^ that small changes in the absolute 
value of the current when the atoms are relaxed are ob- 
served even at large biases. Similar results are valid for 
the other wires as well. 

In conclusion, we have reported first-principles calcu- 
lations of current-induced forces in Si atomic wires. We 



find that these forces are generally strongly nonlinear 
as a function of bias. Since the direct force from the 
bare electrodes is almost linear as a function of bias, 
the nonlinearity originates from the competition between 
the scattering-state and discrete-spectrum force. We 
also find that the average force in the wire is larger 
the shorter the wire, suggesting that atomic wires are 
more difficult to break under current fiow with increas- 
ing length. Finally, current-induced relaxations are found 
to change only slightly the absolute value of the current 
even at large voltages. A similar effect has been pre- 
dicted for molecular wires, indicating it is a general trend 
for nanoscale wires when coherent scattering is the main 
transport mechanism. 
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Similar considerations explain the almost linear depen- 
dence of the force on the single-atom wire (see Fig. pj(d)). 
In this case the extra charge on the atom is almost zero 
at any bias with consequent zero contribution from the 
discrete spectrum, and linear contribution from the con- 
tinuum spectrum as a function of bias. 
It is interesting to note at this point that the sum of the 
continuum- and discrete-state contributions to the force is 
larger for the second atom from the right electrode at all 
biases. A similar result has been inferred by Brandbyge et 
al. [M. Brandbyge, N. Kobayashi, and M. Tsukada, Phys. 
Rev. B 60, 17064 (1999)] from the analysis of the voltage 
drop in gold wires, suggesting that this effect is material 
independent. 

Using a tight-binding approach Todorov at al. have shown 
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that in linear response the largest current-induced force is 
almost constant for wire lengths of three atoms or more. [Q] 
We show here that this trend is also valid for the average 
force at any bias. 



FIG. 3: Extra charge (see text) on the three Si atoms in the 
3-Si wire as a function of bias. 



FIG. 1: Total force as a function of bias in atomic wires con- 
taining (a) two, (b) three, and (c) four Si atoms. The inset 
shows a schematic of one of the wires investigated. The atoms 
are labeled with increasing number from the left electrode (see 
inset), (d) The average force in the Si wires. The left elec- 
trode is positively biased. 



FIG. 2: Different contributions to the total force for a 3-Si 
wire: (a) direct force (see text), (b) scattering-state contribu- 
tion, and (c) contribution from the discrete spectrum. 
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FIG. 4: I-V curve of the 3-Si wire with and without the ef- 
fect of current-induced atomic relaxations. The inset shows 
the unrelaxed (open circles) and relaxed (full circles) atomic 
positions of the Si atoms between the two electrodes (vertical 
thin rectangles) at selected biases. The atoms move opposite 
to the electron flow. 
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